Electronic and Magnetic Reconstructions in Manganite Superlattices 
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We investigate the electronic reconstruction at the interlace between ferromagnetic metallic (FM) 
and antilerromagnetic insulating (AFT) manganites in superlattices using a two-orbital double- 
exchange model including superexchange interactions, Jahn- Teller lattice distortions, and long range 
Coulomb interactions. The magnetic and the transport properties critically depend on the thickness 
ol the AFI layers. We locus on superlattices where the constituent parent manganites have the same 
electron density n — 0.6. The induced ferromagnetic moment in the AFI layers decreases mono- 
tonically with increasing layer width, and the electron-density profile and the magnetic structure 
in the center ol the AFI layer gradually return to the bulk limit. The width ol the AFI layers and 
the charge-transfer profile at the interlaces control the magnitude ol the magnetoresistance and the 
metal-insulator transition ol the FM/AFI superlattices. 



I. INTRODUCTION 

Correlated electron materials often involve the com- 
petition between various ordering tendencies of charge, 
orbital, spin, and lattice degrees of freedom. If fur- 
ther supplemented by weak disorder cluster coexistence, 
percolative transport, and colossal response phenomena 
emerge. It therefore remains a continuing challenge 
to understand the functional properties of transition- 
metal oxides (TMOs}ii^. The physics at the surface of 
TMOs is further enriched by atomic and electronic re- 
constructions and complicated by the lack of inversion 
symmetry^. When surfaces of two TMOs were joined 
together to form an interface^, new phenomena were dis- 
covered in the last decade. In some cases the phases 
at the interface are not even realized in the bulk of ei- 
ther of the TMOs which are joined together— i^. For 
instance the discovery of the two-dimensional electron 
liquid^ which forms at the interface between the insula- 
tors LaAlOa and SrTiOa started a new subfield in the 
research on oxide interfaces. Subsequent experiments 
revealed superconductivity^, ferromagnetisni^, and even 
their coexistence^. These unexpected phases at interfaces 
pose fundamental physics questions and simultaneously 
bear promises of technological importance for the design 
of novel materialsi^. 

The perovskite manganese oxides are a particularly 
remarkable example for the mutual coupling of elec- 
tronic and lattice degrees of freedoroii"— . The man- 
ganites REi_a:AEa;Mn03, where RE and AE denote rare 
and alkaline earth elements, respectively, are known for 
their colossal magnetoresistance^. The various phases 
in manganites with charge, orbital, and magnetic order 
have been elaborated for different combination of RE and 
AE elements and doping regimes ce^iiiii^ii^. The recent 
development in superlattices has created yet another tool 
to explore electronic and magnetic phases of manganites. 

Lai_:rSr^Mn03, a solid solution of LaMnOa (LMO) 
and SrMnOs (SMO) exhibits a number of phases de- 
pending upon the doping concentration x. In addition 
to A-type at cc = and G-type antiferromagnetism at 



X = 1— also ferromagnetic and C-type antiferromag- 
netic (AF) phases exist at low temperature in an inter- 
mediate doping regime. For example Lao.eTSro.saMnOa 
is a ferromagnetic metal and Lao.saSro.eyMnOs is an 
AF insulator— li^ii^. In manganite superlattices 2m lay- 
ers of LMO are deposited on m layers of SMO ^°i^^ 
or vice veiss^. For m = 1, LM02,ri/SM0,„ and 
LM0,„/SM02m are the superlattice counterpart of the 
Lao.erSro.saMnOa and Lao.asSro.eyMnOs manganites, re- 
spectively. 

For LM02m/SM0„i superlattices the uniform FM be- 
havior of the solid solution is recovered for m < 2. 
A metal-insulator transition (MIT) occurs for m > 3 
accompanied by a strongly modulated magnetic struc- 
ture across the interface^i^. Similarly, LM0m/SM02m 
shows a C-type AF phase for m = 1 and m = 2 that exist 
in the solid solution counterpart Lao.saSro.eTMnOa. The 
Neel temperature 250 K of the solid solutions increases to 
320 K for m = 1. For m < 2, charge transfer from LMO 
to SMO may retain properties similar to the solid solu- 
tions. In these two scenarios, the charge transfer is likely 
confined to 1-3 unit cells. For m > 2, the increasing dis- 
tance between the interfaces limits the charge transfer to 
the near vicinity of the interface. Due to this charge con- 
finement LMO and SMO regain their respective parent 
phases away from the interface. 

Considerable progress has been achieved in describing 
the modulated magnetic structure^^"— and the MIT— in 
the LMO/SMO superlattices. The standard two-orbital 
double exchange model for manganites was implemented 
to calculate the average electron density for each layer 
of the LMO/SMO superlattices; away from the interface 
LMO and SMO layers regain their initial electron den- 
sities. The electron-density profile across the interface 
is determined by electrostatics with an electron transfer 
from the LMO to the SMO side. The magnetic-profile 
for each layer follows the bulk phase diagram at the cor- 
responding local density along with a weak proximity 
effect23^. 

In the initial studies of manganite superlattices^^Ti^i 
the primary intent was to improve the tunneling 
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magnctoresistance22iS in a FM/AFI supcrlattice where 
the constituent parent FM and AFI manganites have the 
same electron density, n ~ 1 — x. FM/AFI superlatticcs 
with different combinations of manganites (FM = LSMO, 
LCMO; AFI = PCMO, GdCMO; n = 0.67)^'^-^ have 
been studied in recent years. The observed induced fer- 
romagnetic moment in the AFI layers depends on their 
width and the induced moment is tunable by an external 
magnetic field^'^'^. Remarkably the required external 
field is much smaller than the magnetic field required 
to induce ferromagnetism in the parent AFI compound. 
This results in a large magnetoresistance in the FM/AFI 
superlatticcs which are therefore candidate materials for 
colossal magnetoresistance at room temperature. Man- 
ganite trilayers can also be used to design an efficient 
spin valvG^. 

In this paper, we have studied in detail the electronic 
and magnetic reconstructions in FM/AFI superlattices 
at electron density n = 0.6 for different widths of the 
AFI layers. Electrons are transferred from the the FM 
to the AFI layers at the interface even though the ini- 
tial electron density in the bulk FM and bulk AFI are 
equal. Our study confirms the underlying one-to-one 
correspondence between the density profile and the mag- 
netic profile in the superlattices. With increasing AFI 
layer width, the induced ferromagnetic moment in the 
AFI layers decreases monotonically and drives an MIT 
in the superlattices. In addition to the large magnetore- 
sistance our calculations explain the observed tempera- 
ture crossover between positive and negative magnetore- 
sistance and the two ferromagnetic transitions observed 
in the superlattices^i^. 

The paper is organized as follows: In Sec. II, we intro- 
duce the two-orbital model for bulk perovskite mangan- 
ites. For superlattices we specify the essential modifica- 
tion by adding long-range Coulomb (LRC) interactions 
and briefly present the applied Monte Carlo technique. 
The parameter space of the FM/AFI manganite super- 
lattices is discussed in Sec. III. Electronic and magnetic 
reconstructions at manganite interfaces are emphasized 
in Sec. IV while Sec. V is devoted to the MIT. Resuhs 
with and without LRC interactions are compared in Sec. 
VI. In Sec. VII, various combinations of clcctron-phonon 
couplings for FM and AFI manganites are considered. 
The role of disorder at the interface is discussed in Sec. 
VIII, and conclusions are presented in Sec. IX. 



II. REFERENCE MODEL FOR MANGANITES 

We construct a two-dimensional model Hamiltonian 
for manganite superlattices composed of FM regions sep- 
arated by AFI regions as 

H — HpM + Hafi + Hire, (1) 



where both HpM and Hafi have the same reference 
Hamiltonian^"— 

a/3 

{ij)a i (ij) 

i i iaa 

Here, c and are annihilation and creation operators 
for itinerant eg electrons and a, /? refer to the two Mn- 
Cg orbitals d^ja-ys and d^^^^ri labelled as a and 6, re- 
spectively. The kinetic energy part involves the nearest- 
neighbor hopping of Cg electrons with amplitude 

( ^aa — ^aa = ^: ^bb ^ Hb = ^ab ^ Ka — I V^i ^ab ^ 

= where x and y denote the spatial direc- 

tions on a square lattice. Jh is the Hund's rule coupling 
between t2g spins Si and the Cg electron spin cr,;, and J is 
the AF supercxchange between the t2g spins. A measures 
the strength of the electron-phonon coupling between 
the Jahn- Teller (JT) distortions, Q2i,Q3i and the or- 
bital pseudo spin operators r,2 = J2A4aaCiba + 4ba'^iaa)^ 

= Y.ai4aaCiaa " 4bcr'^iba)- Hcrc Qzi and Qsi are the 
distortions corresponding to the normal vibration modes 
of the oxygen octahedron that remove the degeneracy of 
the Cg levels. The stiffness of the JT modes is denoted 
by K. The stiffness of the breathing mode distortion 
{Qii), which couples to the local electron density is con- 
siderably larger than K and therefore, in the adiabatic 
approximation the coupling to Qu is neglected^ in the 
Hamiltonian Eq. 2. 

Wc treat all t2g spins and lattice degrees of freedom as 
classical^^i^, and measure energies in units of the Mn- 
Mn hopping taa = t. In manganites t is approximately 
0.2-0.5 eV— The estimated value of the Hund's cou- 
pling is 2 eV— , i.e. much larger than t. For this rea- 
son wc further adopt the limit Jh — > oo, for which the 
eg electron spin perfectly aligns along the local t2g spin 
direction. The infinite Hund's rule coupling then natu- 
rally leads to redefine the spinless eg electron operator as 
Cia = cos{9i/2)cia^ + sin(0j/2)e~*'^' Cia|, where the polar 
angle 6i and the azimuthal angle determine the t2g 
spin orientation. In terms of the redefined e^ electron 
operators, the kinetic energy takes the simpler form 

where t^^p is defined as = ^ij^ais with Qij = 
cos(6l,/2)cos(6'j/2) -hsin(6'i/2)sin(6'j/2)e-*('^'-'^^). The 
factor controls the magnitude of the hopping amplitudes 
due to the different orientation of t2g spins at sites i and 

We set if = 1 without loss of generality, because if the 
lattice variable Q is replaced by VkQ, A can be simply 
renormalized as Xj^fK. The length of the t2g spins is 
set to |Si| = 1. In an external magnetic field a Zecman 
coupling Hmag = ^^'^i is addcd to the Hamiltonian. 



The reference 'manganitc model' Href in 2D is con- 
structed to reproduce the correct sequence of magnetic 
phases in the bulk limit^i^. The different phases that 
are captured at low temperature are an orbitally ordered 
insulator aX x = 0^, a charge and orbital ordered ferro- 
magnetic insulator at a; = 0.25^, a ferromagnetic metal 
window around x = 0.67, the CE charge and orbital or- 
dered insulator around x ~ 0.50, a magnetic two dimen- 
sional A-type AF phase, and a G-type AF phase at a; = 
1.00. 

The average electron density of the constituent FM 
and AFI manganites in the FM/AFI superlattices is fixed 
by choosing the same chemical potential throughout 
the superlattice. The LRC interaction between all the 
charges, essential to control the amount of charge trans- 
ferred across the interface is taken into account via a self 
consistent solution of the Coulomb potentials 4)i at the 
mean-field level by setting^ 
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(a) (b) (c) 

FIG. 1: Color online: 1st row: The z components of the 
t2g spins; 2nd row: the electron-density for each site on a 
24x24 lattice at T = 0.01. (a) The parent AF insulator, (b) 
the superlattice with the width w = 11 of the AFI layer, (c) 
the superlattice with w = 11 in the presence of an external 
magnetic field, h = 0.004. 



in the long-range Coulomb part of the Hamiltonian, 



Hire = ^(l^irii. (4) 

i 

It is assumed that all the point charges Zj from the back- 
ground ions are fixed and confined to the Mn sites, {rij) 
refers to the Cg electron density at the Mn site Rj . The 
long-range interactions between the background point 
charges and the Cg electrons are thereby combined with 
the charge- neutrality condition in the superlattice. Al- 
ternatively the Coulomb potentials (pi can be determined 
self consistently by solving the Poisson equation^ which 
is equivalent to solving the Eq. 3. 

The Coulomb interaction strength is controlled by the 
parameter a = /eat where e and a are the dielectric 
constant and the lattice parameter, respectively. This 
mean-field level set up is the minimal basis to study the 
charge transfer across the interfac o^'^'^^i^^ ; it is also com- 
monly used in the context of semiconductor interfaces^. 

The background dielectric constant e is order of 20 in 
manganites^i^. It includes the lattice and the atomic 
polarizability contributions only. The dipolar contribu- 
tion from the charge carrier motion and the associated 
screening is neglected to treat the absolute permittivities 
of FM and AFI manganites on equal footin^^. For a 
lattice constant a — 4A, t = 0.5 eV, and e — 20, the 
approximate value of a is 0.2. The screening length in 
two dimensions (2D) is larger than in three dimensions 
(3D)^. This is taken into account by choosing a larger 
dielectric constant in the 2D model ansatz. Specifically 
we use a = 0.1 in our calculations. The a values em- 
ployed earlier in 3D model calculations were in the range 
0.1 - 1^3^25,26 ^jjjjg much smaller value of a was used 
in ID^. 



We apply an exact diagonalization scheme to the itin- 
erant electron system for each configuration of the back- 
ground classical variables of the i2g spins and the lattice 
distortions. The classical variables are annealed by start- 
ing from a random configuration. At each temperature 
the annealing requires at least 2000 system sweeps, by 
visiting every site of the lattice sequentially and updating 
the system by a metropolis algorithm. For each system 
sweep the computation time scales as 0(A^^) where N is 
the number of lattice sites. With thousands of anneal- 
ing sweeps at various temperatures this procedure is time 
consuming. 

We instead resort to a Monte Carlo sampling technique 
based on the 'traveling cluster approximation' (TCA)^. 
The TCA uses a moving cluster of size Ne constructed 
around the site of the Monte Carlo update. The compu- 
tational cost thereby decreases to 0{NN^) and allows to 
treat system sizes up to ~ 40 x 40, with an 8 x 8 moving 
cluster. The method and the associated transport calcu- 
lations for N = 24 x 24 were extensively benchmarked and 
successfully applied in several earlier studies^i^i^. At 
each system sweep we additionally solve for the Coulomb 
potentials 4>i in Hire self consistently until the electron 
density at each site is converged. The maximum relative 
error for the convergence of the electron density is set to 
0.01 at low temperatures while it is relaxed up to 0.03 at 
higher temperatures. For each fixed set of parameters, 
the calculations start from ten different initial realiza- 
tions of the classical variables. All the physical quanti- 
ties are averaged over the results for these ten starting 
configurations. 
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FIG. 2: Color online: Bulk system at n — 0.6 (a) Distribu- 
tion function for the lattice distortions P(Q) of the bulk API 
(A/ = 1.75) and bulk FM (Aa/ = 1.50) at T = 0.01. (b) The 
resistivity variation with temperature for AFI and FM bulk 
systems. 



III. SUPERLATTICES OF FERROMAGNETIC 
METALS AND ANTIFERROMAGNETIC 
INSULATORS 

Here we analyze specifically superlattices composed of 
FM and AFI manganites. Even with the above described 
simplifications in the model Hamiltonian we have to con- 
tend with manganite states at different doping x, AF su- 
perexchange strengths J, and different electron-phonon 
couplings A. Also the proper choice of relevant parame- 
ters for the combination of ferromagnetic metals and AF 
insulators is not obvious. In addition, due to the slight 
structural and chemical mismatch between the parent 
FM and AFI manganites the parameters at the interface 
may be altered with respect to the bulk values. 

To constrain the parameter space we select the FM and 
AFI manganites of equal electron density n — 1 — x. We 
choose n = 0.6 above half filling to address the existing 
experiments on manganite superlattices near x = 1/3. 
The AF superexchange couplings may vary depending 
on the different combinations of RE and AE ions. But 
for simplicity we will use the typical value J = 0.1— for 
both the FM and the AFI manganites. With the choice 
of the electron density and the superexchange couplings 
we are left with the crucial parameter A to differentiate 
between a FM and an AFI phase. 

Alternatively A can be fixed and the differentiation be- 
tween the ferromagnetic metal and the AF insulator is 
achieved by varying J. This requires a larger value of 
J in the AF insulator, but it leaves the ambiguity for 
choosing modified superexchange couplings at the inter- 
face. Yet another possibility is to keep each parameter 
J and A fixed and to use different hopping amplitudes 
on the two sides of the interface. Again we are left with 
the possible modifications for the hopping amplitude at 
and across the interface. Since the A and J values are 
measured in units of t, a smaller t value is equivalent to 
the combinations of larger A and J values or vice versa. 
Smaller and larger A values are thereby closely related 
to the larger and smaller bandwidth manganites, respec- 
tively. 




CO 

FIG. 3: Color online: Bulk system at n = 0.6 (a) DOS of the 
bulk AFI (A/ = 1.75) and FM (\m = 1.50) at T = 0.01. (b) 
Schematic energy diagram. 

Here we consider two A values to differentiate between 
the FM and the AFI manganites. For the parameters 
,/ = 0.1 and the density n = 0.6, the groundstate is FM 
for A = Xm = 1.50 while it is an AFI for A = A/ = 1.75. 
The AFI phase at n = 0.6 is not a perfect charge and 
orbital ordered CE phase which is stable only at rt = 
^40^42 rpj^g j^Yj p];ja,se at n = 0.6 is inhomogeneous, 
it contains the CE-type zigzag ferromagnetic chains and 
local charge ordered regions as shown in Fig. 1(a). 

We compute the 'one point' distribution function of 
the lattice distortions, P{Q) = {j^ ^{Q~ \Qi\)) where 
N is the total number of sites, to compare them with the 
electron-density variations. The two peaks in Fig. 2(a) for 
the AFI at T = 0.01 indicate that the distortions have 
a bimodal distribution. The sites with a large distortion 
attract more electrons with a density increase to n^ ^ 
0.8. So the spatial density pattern in Fig. 1(a) reflects 
simultaneously the large and small structural distortions. 
The single peak structure of P(Q) for the FM indicates 
instead that there are distortions of similar strength on 
all sites, which facilitates the motion of the electrons and 
thereby supports the metallic character. 

We obtain the resistivity for both FM and AFI bulk 
phases by calculating the d.c. limit of the longitudi- 
nal conductivity as determined by the Kubo-Greenwood 
formulai^ 

q'(^) = ^ ^{n^ - na) ^^'^^ 5{uj ~ {efi ~ e^)), (5) 

a, 13 

with A ~ ire'^ /ha and Ua = 0{fj, — e^), a is the lattice 
spacing, fa/3 denotes the matrix elements of the para- 
magnetic current operator = itj^i ai'^l+x crCi.o- — h.c.) 
between eigenstates \ipa), IV'/s); ^nd ea, are the corre- 
sponding eigenenergies. We extract the d.c. conductivity 
by calculating the 'average' conductivity for a small low 
frequency interval Aoj defined as 

(Tav{Auj) = — / a{u})duj. (6) 
Jo 

Aw is chosen two to three times larger than the mean 
finite-size gap of the system as determined by the ra- 
tio of the bandwidth and the total number of eigenval- 
ues. This procedure has been benchmarkcd in a previous 
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FIG. 4: Schematic view of the FM/AFI superlattice on a 
24x24 lattice. We consider periodic boundary conditions in 
both directions. 



work^. Fig. 2(b) shows the temperature dependence of 
the resistivity p = in units of ha/ne'^. The insulat- 
ing behavior in the AFI at low temperatures is due to 
charge and orbital order, which opens an energy gap in 
the spectrum. 

The density of states (DOS) N{uj) = (^ J2a '^(w-ea)) 
is shown in Fig. 3(a). The center of the gap in the DOS 
of the AFI state is chosen as the energy zero. With this 
choice the Fermi energy of the FM state at the same 
density is at ep = 0.1. The DOS is finite for the FM 
state at its Fermi level. Fig. 3(b) translates the results 
for the DOS into a schematic energy diagram. When the 
FM and the AFI manganites are joined together band 
bending near the interface is expected to shift electrons 
from the FM to the AFI side. An FM/AFI superlattice 
is shown schematically in Fig. 4. The width of the middle 
AFI layer, sandwiched between the FM layers is denoted 
by w. Periodic boundary conditions are enforced in both 
directions and thereby represent a superlattice structure 
composed of alternating FM and AFI regions. 



IV. CHARGE TRANSFER AND MAGNETIC 
RECONSTRUCTION AT THE INTERFACE 

As we show below, if the width of the AF regions is 
small, they may become ferromagnetic and metallic in 
the superlattice structure. The induced magnetization 
is restricted to the near vicinity of the interfaces. The 
distance from the interface at which the parent AF insu- 
lating character is recovered depends upon the parameter 
a and the electron-phonon couplings Xm and A/ . 

In this section we start with the specific choice of 
X]\f = 1.50, A/ ~ 1.75, and a = 0.1, and discuss re- 
sults for the decrease of the magnetization with increas- 
ing width of the AFI layers. This decrease is quantified 
by the average z component of the t2g spins (5*2)7 in the 
AFI layer, where the angular bracket denotes the ther- 
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FIG. 5: Color online: The average z component of t2g spins, 
{Sz)i and the average staggered charge order {CO)i (see text) 
in the AFI layer at T = 0.01. (b) Ferromagnetic structure 
factor S/(0) for different widths of the AFI layer at T = 0.01. 
The average z component is also re-plotted as the dotted line. 



mal average combined with an additional average over 
ten different 'samples'. We also define a measure for the 
local staggered charge order by (CO)/ = J2ieAFi 
gi(7r,7r) ri f^j. |-j^g jg^ygj.. j dcnotcs One of the Nj lattice 
sites in the AFI layer with position r^. In Fig. 5(a), we 
plot {Sz)i and {CO)i for different widths w. For w < 5, 
the averaged {Sz)i is near unity and starts to decrease 
for w > 5, while {CO)i remains small for w < 5 and 
starts to rise for w > 5. Similarly (see Fig. 5(b)) also the 
long- wavelength magnetic structure factor >S'/(0) starts 
to decrease beyond w = 5, where S'/(q) J2ij£AFi 

(Si • Sj) c**! '-'''"""^^. So the decrease in the magnetization 
beyond a critical width Wc is accompanied by emerging 
charge order in the AFI layer. The critical width for 
fixed A/ may change for a different set of parameters Am 
and a. This will be discussed in sections VI and VII. 
The critical width is therefore not directly related to an 
intrinsic length scale of the parent bulk AFI. 

The averaged magnetization necessarily does not re- 
veal any spatial variation transverse to the interface. In 
order to analyze the magnetization profile, we calculate 
the {Sz{x)) for each line of the superlattice with trans- 
verse coordinate x. Fig. 6(a) shows the line averaged 
{Sz{x)) vs. line index for w = 7. {Sz{x)) decreases for 
X = 12, 13, and 14, i.e. in the center lines of the AFI 
layer. So the induced ferromagnetic moment in the AFI 
layer is indeed confined to just the two lines nearest to 
the interface. 

In the bulk for A/ = 1.75 the AFI phase extends to 
electron densities larger than n — 0.6 and changes to 
a ferromagnetic insulating phase at n ~ 0.75*-. It is 
therefore natural to explore the connection between the 
induced magnetization in the AFI layer with the electron 
densities at the interface. As shown in Fig. 6(b), at 
the interface the line averaged electron density {n{x)) of 
the FM layer decreases while {n{x)) for the AFI layer 
increases beyond the initial electron density 0.6. In fact 
{n{x)) at the fully magnetized lines a; = 10 and 16 is 
nearly equal to 0.75. 

The direction of electron transfer from the FM to the 
AFI regions was anticipated already in the discussion in 
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FIG. 6: Color online: Line averaged (a) z component of the t2g 
spins {Sz{x)}, (b) electron density {n{x)), (c) LRC potentials 
{(j>{x)), and (d) lattice distortions {Q{x)) vs. line index for w 
= 7. Open and closed symbols are for the FM and the AFI 
layers, respectively. The temperature is T = 0.01 



FIG. 7: Color online: Line averaged electron density {n{x)) 
for (a) w = 7 with [a = 0.1) and without (a = 0) LRC 
interaction, (b) w = 5, (c) w = 3, (d) w = 11. Open and closed 
symbols are for the FM and the AFI layers, respectively. The 
temperature is T = 0.01 



Section III. Away from the interface the average elec- 
tron density gradually returns to the initial electron den- 
sity n = 0.6. In addition to the overall charge transfer, 
is spatially modulated perpendicular to the inter- 
face. These charge modulations are Friedel-like density 
oscillations; their amplitude decreases with the distance 
from the interface. The spatial density variations are 
nearly symmetric around the central line, x = 13. As 
we have verified a perfect symmetry can be achieved by 
averaging over a larger number of samples. 

For the incommensurate filling n = 0.6, the parent AFI 
is fragmented into small regions with and without charge 
and orbital order. High-density sites (see Fig. 1(a)) are 
always accompanied by large lattice distortions. Simi- 
larly a one-to-one correspondence between the the elec- 
tron densities and the lattice distortions exists in the su- 
perlattice. The line averaged lattice distortions {Q{x)) 
(see Fig. 6(c)) are infact larger in the AFI layer. {Q{x)) 
is modulated similarly to the line averaged electron den- 
sity shown in Fig. 6(b). 

The amount of charge transfer from the FM to the 
AFI layers is controlled by the LRC interactions. The 
line averaged Hartree potential {(j){x)) in Fig. 6(d) is neg- 
ative for a; < 4 and x > 22. The average electron densi- 
ties for those lines is nearly equal to the initial electron 
density n = 0.6. The charge transfer is restricted to a 
few lines near the interface where the averaged poten- 
tials {(f>ix)) are positive and thereby counteracts further 
electron transfer. Since and Ui are the self-consistent 
solutions of the Poisson equation, the increase or the de- 
crease of the averaged electronic density has a one-to-one 



correspondence to the negative or positive curvature of 
(</)(x)). 

For a ~ 0.1, see Fig. 7 (a), the electron densities are 
clearly lower (higher) in the AFI (FM) layer as compared 
to a = 0. The LRC potentials decrease the electron den- 
sity in the center of the AFI layer, but the average elec- 
tron densities for individual lines at a; = 12, x = 13, and 
X = 14 are still enhanced due to the Friedel-like den- 
sity oscillations. The electron density in the central line 
X ~ 13 results from the superposition of the oscillations 
originating from a; 10 and x = 16. The construc- 
tive superposition is more pronounced for the AFI layer 
width w = 5 shown in Fig. 7(b), for which the center line 
is fully magnetized due to the enhanced electron density 
and the induced magnetism originating from the neigh- 
boring magnetized 2+2 interfacial lines. For w = 3, in 
Fig. 7(c), the electron density of the central line is nearly 
equal to the initial electron density and the charge trans- 
fer is truly confined to the interface. For w > 7, density 
oscillations in the AFI layer disappear due the recovery 
of the bulk AFI character (see Fig. 7(d) for w ~ 9). 



V. METAL-INSULATOR TRANSITION 

In this section we compare the magnetic structure fac- 
tor iS'/(0) in the AFI layer with the resistivity for the su- 
perlattice for current fiow perpendicular to the interface. 
Fig. 8(a) shows the temperature dependence of Sj{0) in 
the AFI layer for different layer width w. For all temper- 
atures, the ferromagnetic structure factor and the onset 
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FIG. 8: Color online: Temperature dependence of (a) the 
ferromagnetic structure factor 5/(0) and (b) the resistivity 
p for different widths of the AFI layer. The dotted lines in 
(a) and (b) are for the bulk FM. (c) Si{0) and (d) resistivity 
variation with temperature for w = 11 in the presence of the 
external magnetic fields h = 0.002 and 0.004. Insets in (b) and 
(d) show the resistivity variation in the linear temperature 
scale. 



temperature for ferromagnetism decrease with increasing 
w. The z component of the t2g spins (Szi) and the elec- 
tron densities for each site in the superlattice are shown 
in Fig.l (b) for w = 11. The two lines in the AFI layer 
nearest to the interface are aligned ferromagnetically to 
the t2g spins on the FM side in accordance with the dis- 
cussion in Section IV. The magnetic structure of the cen- 
ter lines is already similar to the bulk AFI phase, and the 
density profile reveals the emergence of local charge or- 
dering patterns in the AFI layer. The superlattice for w 
= 11 is therefore likely to be an insulator. 

Fig. 8(b) displays the temperature dependence of the 
longitudinal resistivity p of the superlattice. p for w = 
11 steeply rises towards low temperatures as expected for 
an insulator. For w < 7, the superlattice is ferromagnetic 
and metallic at low temperatures. For w = 5, the onset 
temperature for ferromagnetism in the AFI layer is near 
0.05, and simultaneously the resistivity start to decrease. 

The low-temperature resistivity increases 
monotonously with the AFI layer width. An MIT 
occurs beyond a threshold width similar as in the 
experiments by Li et al^ who measured the variations 
in the resistivity in LCMO/PCMO supcrlattices for 
different widths of the PCMO layer. In addition, the 
inset in Fig. 8(b) shows a hump in the resistivity for w 
= 9 and w = 7 at the ferromagnetic onset temperatures 
T~ 0.04 and 0.05, respectively, while there is no such 
hump in the bulk FM manganite (see the dashed line of 
Fig. 8(b)). The humps in the resistivity curves increase 



and shift to lower temperatures with increasing AFI 
layer width (for w < 9). This trend agrees well with 
the LCMO/PCMO superlattice experiment!^ with the 
exception of the experimentally observed hump in the 
bulk FM limit. The temperature dependence of p near 
the MIT in bulk FM manganites is tied to the presence 
of intrinsic inhomogencities and disorder— i^. Here we 
have not included any disorder— in the Hamiltonian; the 
MIT in the FM/AFI superlattice has a different origin. 
It is due to the induced onset of ferromagnetism in the 
AFI layer. 

The natural next step is to explore the effect of an 
external magnetic field on the MIT in the superlattice. 
Fig. 8(d) shows that the resistivity for w = 11 is lowered 
by the magnetic field at low temperatures. The connec- 
tion between the resistivity decrease and the onset of fer- 
romagnetism in the AFI layer is verified in Fig. 8(c). The 
MIT for h = 0.004 is at T ~ 0.05; above this temperature, 
the resistivity is larger than the resistivity calculated for 
h = 0. The magnetoresistance, [p{h) — p{0)]/p{0) there- 
fore changes sign near the onset temperature for ferro- 
magnetism. 

At higher temperatures, the t2g spins are randomly 
oriented. We recall that in the limit Jh — > oo the Cg 
electron spins are perfectly aligned along the local t2g 
spin direction. So the current is equally carried by up- 
spin and down-spin electrons. In the presence of a weak 
external magnetic field, the number of eg electrons is in- 
creased in the up-spin channel only in the FM layers and 
the tunneling current from the down-spin channel is de- 
creased. The up-spin channel in the AFI layer remains 
unaltered and restricts the possible enhancement of the 
tunneling current from the up-spin channel. This results 
in an overall increase of the resistivity and a positive 
magnetoresistance. 

At low temperatures in an external magnetic field the 
t2g spins in the AFI layer tend to align in the same di- 
rection as the t2g spins in the FM layers. So the tun- 
neling current in the down-spin channel decreases but 
it increases for the up-spin channel. This decreases the 
overall resistivity of the superlattice. The corresponding 
crossover of the magnetoresistance from negative to pos- 
itive was indeed reported in Lao.rCeo.sMnOs/SrTiOs- 
Nb^ and Lao.82Cao.i8Mn03/SrTio.95Nbo.o503^ hetcro- 
j unctions. 

The z component of the t2g spins and the electron 
densities for each site in the superlattice arc shown in 
Fig. 1(c) for w = 11 and h = 0.004. Ferromagnetic and AF 
structures coexist spatially separated in the AFI layer. 
The ferromagnetic phase in the upper panel is directly 
connected to the charge disordered regions in the lower 
panel. A magnetic field h > 0.007 is required to tune 
the parent AFI into a ferromagnetic metal. Lower fields 
are instead required for the FM/AFI supcrlattices which 
are therefore candidates to design materials with large 
magnetoresistance. 
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FIG. 9: Color online: For a = (a) ^(O) for different widths 
of the AFI layers at T = 0.01. S/(0) for a = 0.1 is also re- 
plotted as the dotted line. Temperature dependence of (b) 
5/(0) and (c) the resistivity for different widths of the AFI 
layers. Legends in (b) and (c) are the same. Inset in (c) 
shows the temperature dependence of the resistivity in the 
linear scale, (d) The temperature dependence of ferromag- 
netic structure factor for the AFI layer, the FM layer for w = 
11, and the bulk FM. 

VI. RESULTS WITHOUT LRC INTERACTIONS 

In order to analyze the effect of the LRC interactions 
we examine the ferromagnetic structure factor S/(0) for 
a = in Fig. 9(a). The critical width, beyond which S/(0) 
starts to decrease, is Wc = 7. With the LRC potentials, 
a = 0.1, the critical width Wc = 5, as discussed in section 
IV, is smaller due to the decrease in the average electron 
density in the AFI layer which is evident in Fig. 7(a). All 
the discussions in the previous two sections remain qual- 
itatively valid also for a = 0, but with a larger critical 
width Wc. 

In Fig. 9(b), we plot the temperature dependence of 
5/(0) in the AFI layer for a = 0. Similar to the results 
in Fig. 8(a), the onset temperature for ferromagnetism 
decreases with increasing w. Fig. 9(c) shows the temper- 
ature dependence of the longitudinal resistivity of the su- 
perlattice. For w = 13 and w = 15, the low-temperature 
resistivity rises by orders of magnitude. The inset of 
Fig. 9(c) zooms into the MIT at intermediate tempera- 
tures on a linear scale. The humps in p(T) decrease and 
shift to higher temperatures for w < 11. 

The onset temperature of ferromagnetism is higher in 
the FM layers than the AFI layers as becomes obvi- 
ous from the temperature dependent structure factor in 
Fig. 9(d). So there are two ferromagnetic transitions in 
the superlattice: the first at Td for ferromagnetism in 
the FM layers and the second for the global fcrromag- 
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FIG. 10: Color online: Si{0) for different combinations of \m 
and \i values at T = 0.01. (b) Si{0) without (A = 0) and 
with (A = 0.2) quenched disorder, and with magnetically 
disordered interfaces (MDIs) (see text) for \m ~ 1.50 and 
Xi — 1.75 at T = 0.01. The dotted line join the three points 
namely w = 1 (for MDIs), w = 5 (for A = 0.2), and w = 9 
(for A = 0). 

netism at Tc2 < Td. The rise in the resistivity p near 
the hump just below Td is apparently due to the on- 
set of ferromagnetism in the FM layers. The rise in p 
is expected for the same reason as for positive magne- 
toresistance discussed in section V. The downturn in p 
results from the onset of global ferromagnetism in the 
super lattices at Tc2- These two transition temperatures 
are also observed in the LCMO/PCMO superlattices^. 
The temperature dependence of the ferromagnetic struc- 
ture factor in the bulk FM state is included as the dotted 
line in Fig. 9(d). The ferromagnetic onset temperature for 
the FM layers in the FM / AFI superlattices is lower than 
for the bulk. 



VII. VARIATION OF Xm AND Aj 

In this section wc compare Si{0) for different combi- 
nations of electron-phonon couplings Aj\/ and A/ . In the 
previous sections Aj\/ = 1.5 and A/ = 1.75 were chosen. 
If A/ is kept fixed and Aj\/ varies between 1.0 and 1.6, 
the groundstate is ferromagnetic and metallic in the bulk 
limit at n = 0.6. The induced magnetization in the AFI 
layer changes little for Am < 1.5 and for different widths 
w (see Fig. 10(a)). For Am = 1.6 however, Si{0) is re- 
duced for w > 5. The difference results from the decrease 
in the electron density in the FM layers nearest to the in- 
terface as becomes evident from the comparison with the 
bulk phase at n = 0.5^. For Am < 1-5, the groundstate 
at n = 0.5 is a ferromagnetic metal while it is a CE-type 
antiferromagnet for Am = 1-6. For Am ^ 1-5, the magne- 
tization in the interfacial FM layers is not altered, even if 
the electron density decreases towards n = 0.5. But for 
A A/ = 1.60, ferromagnetic order becomes unstable at the 
interfacial lines. The smaller magnetization at the inter- 
facial FM layers decreases the induced magnetization in 
the AFI layer for w > 5. 

Fig. 10(a) also indicates the reduced magnetization in 
the AFI layer for the larger electron-phonon coupling 
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\j = 1.8, for which the AFI layer recovers the AF, charge 
ordered state at a smaller width w. It is therefore easier 
to induce a ferromagnetic moment in large-bandwidth 
(small A) manganites. At n = 0.67, the bandwidth of 
PCMO is the largest among those manganites for which 
an AF, charge ordered insulating phase is experimen- 
tally observedi^. In fact, in most of the experimental 
FM/AFI superlattices at n = 0.67 the insulating man- 
ganite PCMO is used along with a variety of different 
FM manganites. 

VIII. DISORDER AT THE INTERFACE 

The size mismatch of RE and AE elements in the man- 
ganites REi_a:AEa:Mn03 leads to tilting and distortions 
of the MnOe octahedra and variations in the local elec- 
tronic parameters. The tilting and distortions of MnOg 
octahedra is generally known as A-type disorder. Here 
we have neglected the intrinsic A-type disorder in both 
the FM and the AFI manganites. Their interface is more 
prone to disorder due to chemical intermixing, lattice 
mismatch, and A-type disorder. We test the effect of 
quenched binary disorder in the terminating line of the 
FM layers at the interface by adding the potential disor- 
der term by ejrij to the Hamiltonian. The sum over 
j is restricted to the terminating lines of the FM layers 
and Cj is the quenched disorder potential with Ij = and 
values ±A. 

In Fig. 10(b), we show the induced magnetization in 
the AFI layer for A = 0.2 and different widths w, which 
apparently changes little between A = and A = 0.2. 
All the discussions in the previous sections remain qual- 
itatively unchanged for moderate disorder A = 0.2. The 
magnetization of the interfacial lines of the FM layers, 
where A is included, is nevertheless likely to decrease for 
larger A. Ultimately the magnetization in those interfa- 
cial lines will be quenched due to spin disorder. In order 
to see the effect of the magnetically disordered interfaces 
on the AFI layer we have fixed randomly oriented spins 
in the interfacial lines throughout the Monte Carlo sim- 
ulations. Indeed, in this extreme case, the induced mag- 
netization in the AFI layer is very small (see Fig. 10(b)), 
irrespective of the AFI layer width. The magnetization 
of the FM layer at the interface is therefore crucial to 
induce the ferromagnetic moments in the AFI layer. 

In the experiments it is apparently not clear whether 
or not there are magnetically disordered interfaces on 
the FM sides in the FM/AFI superlattices^. Thin AFI 
layers are likely to be more susceptible to disorder due 
to strain effects. The disorder strength A at the inter- 
face will therefore increase with decreasing the AFI layer 
width; the width dependence of the disorder strength is 
however hard to quantify. If we combine the results in 
Fig. 10(b) such that disorder strength is maximum for 



smaller AFI layers width and decreases thereafter for 
larger w (see the dotted line) then the magnetic mo- 
ment in the AFI layer behaves non-monotonically which 
agrees qualitatively with the LSMO/PCMO superlattice 
experiments^. But more theoretical work is needed to 
understand this non-monotonic behaviour. 



IX. CONCLUSIONS 

We have investigated the magnetic and electronic prop- 
erties of manganite superlattices at the specific electron 
density n = 0.6 using a two-orbital double-exchange 
model including super-exchange interactions, JT distor- 
tions, and LRC interactions in 2D. At the interface elec- 
trons are transferred from the the FM to the AFI man- 
ganites. Friedel-like density oscillations are observed in 
both the FM and the AFI layers. Due to the charge trans- 
fer and the induced magnetization the smaller width AF 
insulator sandwiched between FM manganites turns into 
a ferromagnetic metal at low temperatures. The induced 
magnetization in the AFI layer decreases beyond a criti- 
cal width and away from the interface the AFI layer grad- 
ually returns to the bulk AFI state. This results in an 
MIT in the FM/AFI superlattice driven by the AFI layer 
width. Only weak external magnetic field is required to 
tune the insulating into a metallic state. The sign of 
the magnetoresistance changes at higher temperatures. 
A remarkable result is that colossal magnetoresistance 
in manganite superlattices is achieved at lower magnetic 
fields as in the bulk phases of the constituent materials. 
For the range of AFI layer widths, at which FM/AFI su- 
perlattices show a metallic behavior at low temperatures, 
a hump in the temperature dependence of the resistivity 
p is observed. The hump in p is due to the presence of 
two ferromagnetic transition temperatures in the super- 
lattice. The height of the hump decreases and shifts to 
higher temperatures with decreasing AFI layer width. 

In summary, the width of the AFI layers controls the 
magnitude of the magnetoresistance and the MIT in the 
FM/AFI superlattices. Our 2D model Hamiltonian cal- 
culations provide a basis for explaining the MIT in man- 
ganite superlattices. The non-monotonic behaviour of 
the induced ferromagnetic moment in the AFI layer is 
traced to the effects of magnetic disorder at the inter- 
face. 
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